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We study numerically and analytically a model of self-propelled polar disks on a substrate in two 
dimensions. The particles interact via isotropic repulsive forces and are subject to rotational noise, 
but there is no aligning interaction. As a result, the system does not exhibit an ordered state. The 
isotropic fluid phase separates well below close packing and exhibits the large number fluctuations 
and clustering found ubiquitously in active systems. Our work shows that this behavior is a generic 
property of systems that are driven out of equilibrium locally, as for instance by self propulsion. 

PACS numbers: 



Collections of self-propelled (SP) particles have been 
studied extensively as the simplest model for "active ma- 
terials", a novel class of nonequilibrium systems com- 
posed of interacting units that individually consume 
energy and collectively generate motion or mechanical 
stresses [Ij. Active systems span an enormous range of 
length scales, from the cell cytoskeleton |2j, to bacte- 
rial colonies [3 , tissues [4 and animal groups [5]. These 
disparate systems exhibit common mesoscopic to large- 
scale phenomena, including swarming, non-equilibrium 
disorder-order transitions, pattern formation, anomalous 
fluctuations and surprising mechanical properties [6[ |7] . 

Active particles are generally elongated and can or- 
der in states with either polar or apolar (nematic) ori- 
ent at ional order. A remarkable property of such ordered 
states are giant number fluctuations. In equilibrium sys- 
tems, away from continuous phase transitions, the stan- 
dard deviation in the mean number of particles N 
scales as \fN for ^ oo. In active systems AA^ can 
become very large and scale as A^^, with a an exponent 
predicted to be as large as 1 in two dimensions [7H9] • This 
theoretical prediction has been demonstrated experimen- 
tally [TOHIl] and verified in simulations of agent-based 
models [TSHIS]. Both nematic and polar states exhibit 
giant number fluctuations, which are believed to be as- 
sociated with the broken orientational symmetry. 

In this paper we study a model of SP soft repulsive 
disks with no alignment rule. Since the particles are 
disks, steric effects, although included in the model, do 
not yield large scale alignment, in contrast to SP rods 
that can order in nematic states [161 |T7]. As a result, our 
particles, although self-propelled, do not order in a mov- 
ing state at any density. Figure [l] shows, however, that 
above a packing fraction 0c ^ 0.4 this minimal system 
phase separates into a solid-like and a gas phase, hence 
exhibits giant number fluctuations for (j) > (j)c (Pig- [2]). 
While the giant fluctuations seen in the ordered state 
of nematic and polar active systems [TOl [131 [El are be- 
lieve to be intimately related to the broken orientational 
symmetry, the ones seen here arise in the absence of any 
broken symmetry when the rate at which self-propulsion 
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FIG. 1: (color online) Snapshots of Nt = 10^ disks for = 
0.39 (top row) and = 0.7 (bottom row). Same-size clusters, 
defined by particles overlap, are highlighted by color coding. 
The left frames are for a thermal system at ksT — 0.1. The 
right frames are for SP disks with vq — 1 and i/^ = 5 x 10~^. 



is suppressed due to steric trapping exceeds the rate of 
density convection, resulting in phase separation. Simi- 
lar clustering has been seen in spherical vibrated granu- 
lar particles, although there inelasticity of the interaction 
may play a role [181 [E] • Our work supports the sugges- 
tion by Gates and collaborators [2QH22j that clustering 
and phase separation are generic properties of systems 
that are driven out of equilibrium by a persistent local 
energy input that breaks detailed balance. We further 
show that the notion of effective temperature cannot be 
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FIG. 2: (color online) Standard deviation AN versus the 
average number N of particles in a subsystem of size i = 
y^Tia'^N/cf) for packing fractions (f) from 0.1 to 0.8, Nt = 10^ 
and 1 < N < Nt- The dashed and dotted lines correspond 
to AN — N'^^'^ and AN — N ^ respectively. Inset: exponent 
a such that AN - N"" versus (j) for Nt from 200 to 10^. 



used to describe the system. Finally, using input from 
the numerics, we construct a continuum model that re- 
produces the results of the simulations. 

Numerical Model We restrict ourselves to systems 
without momentum conservation, such as granular mate- 
rials or living organisms on a substrate. The particles are 
soft repulsive disks of radius a with a polarity defined by 
an axis = (cos^^, sin^^), where i labels the particles. 
The dynamics is governed by the equations 



+ r,J{t), dtei = rji{t), (1) 



with vq the self-propulsion speed and /i the mobility. 
The translational and rotational noise terms, rjf{t) and 
77i(t), are Gaussian and white, with zero mean and cor- 
relations {vla{t)vji3{t')) = 2DSijSa(3S{t - t') (the Greek 
labels denote Cartesian coordinates) and {r]i{t)r]j{t') = 
2urSijS{t — t')^ with D = kBTji the Brownian diffu- 
sion coefficient and Vr the rotational diffusion rate. For 
Brownian particles of size a, D ~ 0?^^ at low density. 
Here, however, we treat D and Vr as independent noise 
strengths, with D controlled by thermal noise and Vr a 
measure of nonequilibrium angular noise as it may arise 
for instance from tumble dynamics of swimming organ- 
isms [23^. In the numerical work described below we ne- 
glect the translational noise to highlight the crucial role 
of the angular noise. The force F^j between disks i and 



j is short-ranged and repulsive: 
if Vij < 2a and F^j 



-k(2a 



otherwise. Although the self- 
propulsion speed ^0 is fixed, the instantaneous speed 
of each particle is determined by the forces due to the 
neighbors, in contrast to Vicsek-type models (but see 
Refs. [24H27 1). We have performed molecular dynam- 
ics simulations of Eqs. ([T]) at T = with Nt = 100 to 
10000 particles and periodic boundary conditions in a 



box of size L. We explore the phase diagram by vary- 
ing the self propulsion speed and the packing fraction 
(j) = Ntt^o? / L'^ . In the numerics we have scaled lengths 
with the radius a of the disks and times with 10/ {fik)^ and 
have fixed the rotational noise strength Ur = 5 x 10~^. 

Athermal phase separation. A hallmark property of 
SP particles is strong clustering in the absence of at- 
tractive interactions. This is ubiquitous in the ordered 
state of Vicsek-type models and has been seen even in 
the absence of polar aligning rules in simulations of SP 
hard rods |28H3Q] . In this case the anisotropic shape of 
the particles provides an aligning, although apolar, in- 
teraction, that enhances cluster formation. Our model, 
in contrast, consists of radially symmetric particles and 
no alignment can arise from steric effects. Nonetheless, 
we observe strong athermal clustering due solely to the 
non equilibrium nature of the SP disks. Above a criti- 
cal packing fraction (j)c ~ 0.4, the system separates into 
dense macroscopic clusters and a low density phase. Even 
below (j)c clustering is much more pronounced than in a 
thermal system, as illustrated on Fig. [l] for packing frac- 
tions below (top) and above (bottom) <pc- The frames 
to the right are snapshots of the SP disks model. The 
left frames show snapshots of an "equivalent thermal sys- 
tem", defined as one with comparable overlap between 
particles, at the same packing fractions. In dimensionless 
units, the typical overlap S in the active case is obtained 
by balancing the active force vo//j. with the repulsive force 
kS^ with (5 = 0.1 for ^0 = 1- The equivalent thermal sys- 
tem is then obtained by setting ksT = kS'^ =0.1. This 
comparison indicates that the SP system cannot be sim- 
ply described by an effective temperature, in agreement 
with [H], and in contrast to what has been argued by 
some authors [SJ, 33 . The inadequacy of the notion 
of effective temperature is also supported by an analysis 
of the cluster size probability distribution (not shown). 
The notion of effective temperature may at best hold in 
very dilute systems. To see this we formally integrate the 
angular dynamics and rewrite Eqs. ([T]) solely in terms of 
translational dynamics as 



dfTi 



(2) 



where the noise ^i{t) has zero mean and variance 



{i^am,0{t'))o,=2[DS{t-t')- 



f e-'^^l*-*'l],5«^% (3) 



where {...)eo includes an average over the initial values 
of the angles. This shows that Eqs. ([T]) are equivalent to 
those for interacting soft disks with non-Markovian noise 
with memory time This noise can be approximated 
as white, with {^ia{t)^i^{t'))eo = 2^kB (T ^ T^) 5(t - 

2 

t')Sai3^ij and ksTe = only for \t - t'\ » 1/ur. 

The effective temperature description will be adequate 
only in the very dilute gas phase, when is shorter 
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than the mean free time. In dimensionless units this re- 
quires (j) < itv'^/vq. The parameter values used in Fig. [l] 
give ttv^/vq 10~^ and /c^Tg = 100, i.e. essentiahy zero 
density and infinite temperature. 

Giant Number Fluctuations and Orientational Corre- 
lations. In the phase separated regime for (j) > (j)c 
observe giant number fluctuations: as shown in Fig. |2] 
the variance AA^ of the fluctuations in the number of 
particles in a subregion of size containing a total 
mean number of particles A^, scales as AA^ ~ A^^, with 
a = 0.95 ± 0.05. This exponent is consistent with value 
a = 1 expected for phase separation in 2d. Orientational 
correlations decay exponentially in both phases and the 
large cluster is stationary, confirming the absence of large 
scale orientational order. Residual correlations exist at 
the surface of large clusters, where particles tend to point 
inward. Particles mainly enter or leave the cluster indi- 
vidually or in pairs. These observations demonstrate that 
orientational correlations, to be discussed in a future pub- 
lication, do not play the central role in controlling phase 
separation. 

Mean-square displacement. To further characterize 
the dynamics, we have evaluated the mean square dis- 
placement (MSD) of a tagged disk, shown in Fig. [s] 
An individual SP particle described by Eqs. Q with 
Yij = performs a persistent random walk (PRW), with 



([Ar(t)]2)=4Z)o 



t- 



(e- 



1) 



, with 1^0 = 



2u, 



and 



a crossover from ballistic behavior with ([Ar(t)]^) ~ 'Uq^^ 
for t « to diffusive behavior, with ([Ar(t)]^) 
4:Dot at long times [36 The PRW form fits the data at 
vanishingly small packing fraction. At non-zero density 
the MSD displacement can still be fitted by a PRW form, 
([Ar(t)n = 4De [t+^(e-^^*-l)], with = 
and Ve density-dependent fitting parameters. The 
bottom inset of Fig. [3] shows the effective self-propulsion 
speed Ve{(j)) obtained from the fits as a function of den- 
sity. For (j) < (j)c^ '^e(</^) is well fitted by a linear form 
'Ue(^) = '^0(1 — with A ~ 0.9 independent of vq for 
the three simulated values {vq = 0.5,1 and 2). In con- 
trast, ^ Ur depends weakly on density. As a result, 
^ I^o(l — as shown in the top inset of Fig. [s] 
Above (j)c^ the MSD is slower than ballistic at short times. 
In spite of this, the PRW form with a linear decrease of 
Ve with (j) still fits surprisingly well. This may be due to 
the fact that above (j)c a large fraction of disks belongs to 
stationary clusters and the MSD is controlled by parti- 
cles in the low density regions. All densities shown here 
are below the crystallization density ~ 0.91, above 
which a bounded MSD is expected [31^ . As pointed out in 
Ref. [21], De{(j)) represents a collective diffusivity renor- 
malized by interactions. This quantity was calculated in 
pO^ by a statistical analysis of run and tumble dynamics 
in one dimension, where an exponential decay of with 
density was obtained. The linear decrease of Ve with den- 
sity found in our model is consistent with the predictions 




FIG. 3: (color online) Mean square displacement of a tagged 
disk versus time for various packing fractions 0, showing 
crossover from ballistic to diffusive behavior. Insets: effec- 
tive diffusion constant De{(j))/Do and self-propulsion speed 
Ve{(l))/vo obtained from the fit to the data as functions of 
for bare SP velocity vq = 0.5,1,2 (diamonds, circles and 
squares respectively). 



of Ref. [20^ for weak interaction strength. The resulting 
suppression of diffusion is the mechanism responsible for 
phase separation and cluster formation. 

Continuum model. We now use the findings from our 
simulations to construct an empirical continuum model 
that captures the dynamics of the system. The anal- 
ysis of the MSD indicates that one of the effects of 
steric repulsion is the replacement of vq in Eqs. ([T]) by 
Ve{(l)). After this replacement, we use standard meth- 
ods [34[ |35] to coarse-grain the microscopic dynamics 
and derive continuum equations for the conserved den- 
sity p(r, t) of active particles and the polarization density 
p(r, t) = p(r, t)P(r, t), with P the orientational order pa- 
rameter. Although our system does not order, the noisy 
angular dynamics of Eqs. ([T]), described in the continuum 
by the coupling to P, is crucial in controlling the behavior 
of the system. A minimal version of the hydrodynamics 
of overdamped SP particles [8 that neglects all convec- 
tive nonlinearities, but is adequate for our purpose, is 
given by \37\ 



dtp 



-V-(^eP-^pVp + fp) 



dtp = -VrV - V(Vep) 



(4a) 



(4b) 



where fp and fp represent Gaussian white noise with zero 
mean and correlations {fpi{r^t)fpj{r'^t')) = 2ApSijS{r — 
v')5{t-t') and {U,{v,t)Uj{v',t'))=2/\5,j6{v-v')6(t- 
t'). The density equation is simply a convection-diffusion 
equation, with advection velocity given by the local self- 
propulsion speed. We include a finite value for the den- 
sity diffusion Dp because even in the absence of trans- 
lational noise in the microscopic dynamics, density dif- 
fusion (and polarization diffusion K) would be induced 
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in the system through interactions. The polarization de- 
cays at rate Vr and is convected by pressure-hke gradi- 
ents ~ ^{vep). The only homogeneous stationary state 
and (4b) is the isotropic state 
= 0. 



described by Eqs. (|4a 
with p 



po and p = 0. To examine the stability 
of this state we consider the dynamics of fluctuations 
Sp = p — po and ^p. Introducing Fourier amplitudes 



(<^/>q,a;,©q,a; 



7.< 



(^p(r,t),V-p(r,t)), 



the linearized equations take the form 



wq'^Sp^ 



q,a; 



(5a) 



(5b) 



< 



po 



where w = Ve{po) + Po<(Po). Since <(po) = (^) 

0, w can change sign, signaling self-trapping. The disper- 
sion relations of the linear modes are easily calculated. 
At small wave vector the dynamics is controlled by a dif- 
fusive mode uj-{q) iVq^^ with V = Dp Vew/vr an 
effective diffusivity. If il' > 0, convective currents associ- 
ated with self-propulsion exceed self-trapping responsible 
for the decrease of Ve- Then P > and the isotropic state 
is stable. When w <{) the effective diffusivity V becomes 
negative for '^el'^l > ^r^p, signaling unstable growth of 
density fluctuations and phase separation. Using the lin- 
ear fit for Ve{(l)) from the numerics, we estimate that w 
changes sign at (/>* = 1/(2 A). If we neglect thermal dif- 
fusion {Dp = 0) the isotropic state is unstable for all 
packing fractions > 0* ~ 0.45, that we identify with 
^ 0.4 found in the numerics. A finite value of Dp 
shifts the instability boundary to higher density. 

Static structure factor. The continuum model can 
also be used to evaluate the static structure factor 6'(q) = 
^(Jpq(t)(5p_q(— t)), which is a direct measure of the spa- 
tial correlations of density fluctuations, with S{q ^ 0) = 

We evaluate S{q) in the region > by comput- 
ing the dynamical structure factor 6'(q, cj) = ;^(|^Pq,a;p) 
from the linearized equations for the fluctuations with 
noise. Then ^(q) = f^6'(q, cj). The noise in the 
density equation does not contribute at small wave vec- 
tors, and we obtain 



Siq) = 



{Dp + K)q^][urV + DpKq^] 



(6) 



with S{0) = ^^2^ • As the instability is approached from 
below V ^ and 5'(0) diverges, with a behavior that 
reminds of that at an equilibrium critical point. The 
decay of correlations is characterized by a length scale 
^ = ^/KDp/{iyrV) - (0 - 0c)"^/^ that diverges at the 
instability. The static structure factor obtained from sim- 
ulations is shown in Fig. [IJ For (j) > (j)c^ S{q) diverges at 
small q and is reasonably well described by S{q) ~ q~^ 
with a ~ 2, consistent with what is expected for a phase 
separated system. Below the transition (0 < 0c) , we 
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FIG. 4: (color online) Static structure factor S{q) for (j) from 
0.1 to 0.8 and Nt = 2000. The dashed curve is 5' - g"^ 
Inset: Rescaled structure factor S{q) for (j) < (j)c (see text). 
The dashed curve is = 1/(1 + ^^). 



fit S{q) at low g to a Lorentzian S{q) = S'o/(l + Q^i^)- 
The inset of Fig. [4] shows a good collapse of the rescaled 
structure factor S{q) = S = S/Sq, with q = q(_ and 
and ^ obtained from the fits. The growth of the correla- 
tion length ^ is not inconsistent with a divergence at the 
transition, but a detailed study of the transition region 
is needed to determine scaling exponents. 

In summary, we have shown that self-propelled par- 
ticles with no alignment exhibit an athermal clustering 
instability to a phase-separated regime well below close 
packing. Above the instability, the system exhibits large 
density fluctuations ubiquitous in active systems. 
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